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A real-time formalism is proposed to incorporate finite temperatures into confined quark systems at the 
example of the Generalised Nambu-Jona-Lasinio model for QCD. This approach allows one to study vari- 
ous properties of the system at T > 0, such as chiral symmetry breaking and restoration, properties of the 
bound-state spectrum, and so on. 
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1. INTRODUCTION 



Oh While the property of colour confinement is expected 
from general considerations and is supported by lat- 
itice calculations, no effective theory exists which de- 
rives confinement directly from the QCD Lagrangian. 
Therefore, in the absence of controllable analytic solu- 
tions of QCD, one has to resort to models to describe 
hadronic phenomenology. However, these models must 
^satisfy a very stringent set of minimal requirements. 
They must be: i) relativistic, ii) chirally symmetric, and 
iii) able to provide a mechanism of spontaneous break- 
'ing of chiral symmetry, iv) they should contain confine- 
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ment, in order to be able to address the issue of excited 
states. To satisfy all the above requirements with one 
. . and the same model, turns out to be a highly nontrivial 
rS 'matter. There is a model, however, which does incor- 
"porate the aforementioned set of requirements. This 
'is the Generalised Nambu-Jona-Lasinio (GNJL) model 
with the instantaneous vector confining kernel [U [5] 
and confinement of quarks guaranteed due to an in- 
stantaneous infinitely rising (for example, linear) po- 
tential, to which extra ingredients can be added like, 
for instance, the colour Coulomb potential. Then chi- 
ral symmetry breaking can be described by standard 
summation of valence quark self-interaction loops (the 
mass-gap equation), while hadrons are obtained from 
the Bethe-Salpeter equation for the multiquark bound 
states. Besides, GNJL is known to fulfill the well-known 
low-energy theorems of Gell-Mann, Oakes, and Renner, 
Goldberger and Treiman, Adler selfconsistency zero, the 
Weinberg theorem, and so on — see Ref. [3] . For highly 
excited hadrons this model predicts effective restoration 
of chiral symmetry [4 in accordance with general expec- 
tations [5]. At this stage, we would like to emphasise 



the universal nature of the above low-energy theorems 
as well as of the chiral restoration in excited hadrons 
regardless of the actual giuonic interactions that would 
eventually be responsible for chiral symmetry breaking. 
Thus, the GNJL model is expected to be quite useful 
to describe all those phenomena which are essentially 
driven by the same global chiral symmetry. 

Incorporation of finite chemical potentials and finite 
temperatures stand as important forward steps in the 
development of this class of models which not only help 
us into solving various puzzles posed by experimental 
results but, as importantly, they may be of assistance 
in predicting new phenomena for experimental testing. 
A number of attempts has been undertaken to include 
finite quark densities into consideration [6l [7] . By con- 
trast, the incorporation of finite temperature appears 
to be a slightly more involved problem. Indeed, one 
may encounter technical problems applying techniques 
such as the Matsubara finite-temperature approach to 
the models, where only real-time calculations are possi- 
ble. Thus we adhere to a different approach, based on 
the formalism of Bogoliubov-Valatin transformations, 
which is a natural mathematical tool to deal with GNJL 
systems. Such an approach to chiral symmetry break- 
ing at r = was successfully used in a series of papers 
[2] and then extended to redefine the theory entirely 
in terms of the quark-antiquark bound states [8]. The 
underlying idea of the method is to consider the true 
thermal vacuum of the system, as a collective (coherent- 
like) phenomenon of zero-temperature chirally symmet- 
ric theory eigenstates. To this end it is necessary to 
build a pseudounitary, temperature-dependent, trans- 
formation relating the thermal vacuum with the zero- 
temperature unbroken vacuum (such a real-time formal- 
ism for not confined systems was suggested in Ref. [9]). 
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For quarks, this transformation leads to the formation 
of new effective objects — dressed quarks which, when 
seen from the zero-temperature chirally symmetric Fock 
space, look like original bare quarks "escorted by" an 
accompanying cloud of quark-antiquark pairs. The cor- 
nerstone of this method, at T = 0, amounts to finding, 
through a variational calculus, the minimal vacuum en- 
ergy of the theory evaluated as a functional of the order 
parameter, called the chiral angle. Then in order to in- 
clude finite temperature effects we have to consider the 
free energy Fvac = £^vac — TS, with S being the entropy 
of the vacuum, instead of the usual E'vac = (0|-ff |0). For 
a given temperature T this minimisation process ensures 
a proper balance between the vacuum energy and the 
entropy of the vacuum, so that the resulting vacuum 
state is stable, with Fvac < ^vac, where F^ac is the free 
energy of the trivial, unbroken, vacuum. 



2. GNJL MODEL AT T = 

We start with a short introduction to the GNJL chi- 
ral quark model [H [2] given by the Hamiltonian 

^ '' d^'xd^y J^ix,t)K^,ix-y)J^{y,t), (1) 

where J^{x,t) = 'il){x,t)^^^'ip{x,t). Having in mind 
that the qualitative results and conclusions do not 
depend on any particular form of the quark kernel 
K^u{x — y), we only require it to be confining and to 
introduce a natural scale. For phenomenological appli- 
cations this scale is to be fixed of order of 300 ~ 400 
MeV. No extra constraints are imposed on the kernel. 

Among possible quark-quark interactions described 
by the Hamiltonian ([T]) we have quark selfinteractions. 
It turns out that these self-interactions are removed by 
the use of an appropriate Bogoliubov-Valatin transfor- 
mation from bare to dressed quarks, which can be con- 
veniently parametrised by means of the so-called chiral 
angle fp {p being the relative momentum in the dressing 
quark-antiquark pairs) [1] [2] : 
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^1 + sin(/3p + {a.p)^\ - sixiifp uo{p), 



y/1 + sinipp - {ctp)y/l - sinipp vq{-p) 



venient to define the chiral angle varying in the range 
— 7r/2 < ipp ^ 7r/2 with the boundary conditions 

= 0) = 7r/2, ^p{p ^ (X3) -> 0. (2) 

The normal ordered Hamiltonian ([1]) becomes: 

H = £'vac+ : H2 : + H4 :, (3) 

and the usual procedure is to demand the quadratic part 
: H2 : to be diagonal or, equivalently, that the vacuum 
energy E'vac should become a minimum. Then the cor- 
responding mass-gap equation, 

6E,,,[ip]/dipp^0, Fvac[<^] = (0M|i?|0M), (4) 

ensures the anomalous Bogoliubov terms b'^d^ and db to 
be absent in : iJ2 :• As soon as the mass-gap equation 
is solved and a nontrivial chiral angle (pp is found, the 
Hamiltonian ^ takes a diagonal form. 



H — Eyac + 



E 

p.a,s 



(5) 



where Ep is the dressed-quark dispersive law and the 
quark creation and annihilation operators 6p^, dpi, 
and dpg are consistent with the Fock space T built on 
top of the nontrivial BCS vacuum |0) = |0[(^]). The 
contribution of the : iJ4 : part is suppressed as \/\/Wc, 
so this completes the diagonalisation of the Hamiltonian 
^ in the quark sector (at the BCS level). 

Since the definition of a one-particle state has 
changed dramatically, the vacuum |0) annihilated by the 
dressed single-particle operators is also different from 
the trivial vacuum |0)o = |0[</3 = 0]). Indeed, as it 
always happens after a Bogoliubov-Valatin transforma- 
tion, the true vacuum, with the minimal vacuum energy, 
contains an infinite set of strongly correlated ^Po quark- 
antiquark pairs [2J, 
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(6) 



where a is the colour index of Nc colours. It is con- 



where = bpl[{cTp)icF2\ss' dpi, with cr's being the 2x2 
Pauli matrices. The operator Cj, creates a ^Pq quark- 
antiquark pair with zero total momentum and the rel- 
ative three-dimensional momentum 2p. The chiral an- 
gle ipp is a solution to the mass-gap equation ([4]) and 
it "measures" the weight of the pairs with the given 
relative momentum, so that the operator exp[(5o — Qo\ 
creates a cloud of correlated pairs, and the BCS vacuum 
|0) has the form of a coherent-like state when seen from 
the point of view of the naive Fock space (that is, 
the Fock space built on top of the trivial vacuum |0)o). 
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The quark Fock space is built over |0) by repeated 
application of quark/antiquark creation operators. 

Using the commutation relations for the quark op- 
erators b and d, one can easily find the following repre- 
sentation for the new vacuum 2 : 



io[^])=n 



(7) 
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(8) 

clearly obey the condition WQp -t- wip + W2p = 1 and thus 
they admit the following natural interpretation: they 
represent the corresponding probabilities of having in 
the new vacuum zero, one, and two quark-antiquark 
pairs with a given relative momentum. Notice that pow- 
ers of the operator higher than two cannot appear 
because of the Fermi statistics for quarks. Then it is 
easy to show that the new vacuum is normalised, 



(0[<^]|0[V3]) = JJ(wOp + Wip + W2p) = 1, 



(9) 



and that it is orthogonal to the trivial vacuum in the 
limit of infinite volume of the space V: 



(0[(^]|0)o =exp 



0. 



(10) 



3. GNJL MODEL AT FINITE TEMPERATURES 

As it was mentioned before, consideration of finite 
temperatures amounts to evaluation of both the vacuum 
energy and the entropy of the system simultaneously 
and then to the minimisation of the free energy with re- 
spect to the order parameter. With the help of Eq. ^ 
one can easily find for the entropy of the vacuum: 



(11) 



p n— 



where the probabihties ujfip are defined in Eq. ([5]). 

Although we shall be working only with the c- 
number (jlip , it is instructive to see how it appears as a 
result of averaging over the vacuum of a local operator, 
iS* = (0[iy9]|-fC|0[93]), which should: i) be built entirely in 
terms of the operators Cp and Cj, — the only build- 
ing blocks with the quantum numbers of the vacuum at 
our disposal; ii) commute with the Hamiltonian ([5|) ; iii) 
comply with the relation (fTTjl . The final result reads: 
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As it was mentioned before, the thermal mass-gap 
equation guarantees the balance between the energy and 
the entropy of the vacuum and can be derived by min- 
imising the free energy operator H — TK averaged over 
the vacuum |0[(^]). Then, with the help of Eqs. ([4]) and 
(fTTj) . one can find for the free energy: 

i^vac[(^] - (0[(y3]|(if - T K)\Q[^]) = E,,c " TS. (12) 

Then the mass-gap equation is 



— , SFvac[^]/Sipp = 6Evs,c[v]/Sipp - TSS[ip]/Sipp = 0, (13) 

which is nothing but a generalisation of the zero-tem- 
perature mass-gap equation (j4|) to finite temperatures. 
It is easy to find by a straightforward calculation that 



SS/6ipp 



sm ifp 



V 2 



cos (fp \n2 . (14) 



Mass-gap equation (fT5|) with the entropy and its 
variation given by Eqs. pT|) and ([U]) is the central result 
of this paper. Its physical interpretation is straightfor- 
ward: '^Pq quark-antiquark pairs are condensed in the 
vacuum to lower the vacuum energy and to break chi- 
ral symmetry. With temperature, some pairs are "re- 
moved" from the vacuum so to absorb the heat that 
has been added to the system. The chiral condensate 
decreases (evaporates) accordingly and chiral symmetry 
gets less broken. The mass-gap equation guarantees the 
proper balance between these two opposite processes. 
Concluding this discussion, let us remind the readers 
that the chiral angle is nothing but the wave function of 
the Goldstone boson (^f^ — sin ipp [HE]) responsible for 
the chiral symmetry breaking. With the proper solution 
of the thermal mass- gap equation in hands, one is able 
to study microscopically the process of the Goldstone 
boson "melting" with the temperature increase. 

The extra terms in the mass-gap equation which 
stem from the temperature are consistent with the 
boundary conditions ([2]) imposed on the chiral angle. 
Indeed, as clearly seen from Eq. 

(^^/'5^p)|^,=./2 ^ iSS/5vp\^^^o = 0- (15) 

Since equation SS/Sipp = does not have other solu- 
tions, the entropy increases monotonously from its min- 
imum at ipp — to its maximum at ipp — 7r/2. Then 
Eq. (|13p embodies two natural limits: in the infinite 
quark mass limit one must have ipp = tt/2 for all p's, 
and indeed we can immediately see that the tempera- 
ture dependence of the mass-gap equation vanishes for 
any finite T and we are left with the zero-temperature 
mass-gap equation. On the other hand, for very high 
temperatures and finite bare quark masses, there is only 
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one solution, that is tpp — 0, which is the chiral restora- 
tion scenario. 

An attractive feature of the GNJL model is that the 
qualitative predictions of the model are robust against 
variations of the quark kernel. We therefore consider, 
as a paradigmatic example, the simplest quark kernel 
compatible with the requirements of confinement, 

Kf,^{x - y) = 3^05^0 V^)(^c - y), Vo{r) = K^r^, (16) 

and study the properties of the thermal mass-gap equa- 
tion (|13p numerically. The harmonic oscillator confin- 
ing potential is known to lead to the simplest nontrivial 
mass-gap equation. Indeed, since the Fourier transform 
of the quadratic potential (fTB|) is given by the Lapla- 
cian of the 3D delta-function, the mass-gap equation is 
differential Switching from the harmonic oscillator 
potential to another form of the potential amounts to 
simple modifications of the quantitative results, while 
all qualitative conclusions remain intact (see, for exam- 
ple, Ref. [ini for a detailed analysis). 

Then the vacuum energy is given by the functional 



Ap = m 



(27r)3 



Ap sin ipp 



Bp cos ifip 



Br, 



1 



d^k 
(2^ 
d^k 
(2^ 



Vo(p- fc)sin^fc, 



(17) 
(18) 



{pk)Vo{p ~ k)coinpk- (19) 



Here Cp = \{Nc ~l/Nc) is the SU{Nc)c Casimir op- 
erator in the fundamental representation and the degen- 
eracy factor g counts the number of independent quark 
degrees of freedom, g = (2s + \)NcNf = 12. Also the 
substitution ^ ^ / (^J'yi was made for the sake of 
convenience. The mass-gap equation takes the form: 

[Ap cos ipp - Bp sin ifp] - (r/2) {5S/5ifp) = 0, (20) 

with 6S/Sipp given in Eq. (|14p or, in the explicit form, 

2<y3p sm2ipp 



^0 



P 

T 



2p sin ipp 

,2 fp' 



2m cos fp 



2 In (ta,n^ + cos tpp In 2 



(21) 



In what follows we consider the chiral limit and set the 
quark current mass m — 0. 

In Fig. [1] we plot solutions of Eq. (PTjl for various 
values of the temperature. From this figure one can see 
that, as was anticipated before, with the increase of the 
temperature, the chiral angle approaches chirally sym- 
metric trivial solution. In addition, in Fig. [21 we plot 
the temperature dependence of the chiral condensate 



S{T) = {qq) = 



dppHm<fp{T), (22) 




Solutions to the mass-gap equation (|21|l for the temper- 
atures r = (the upper curve), 0.2Ko, OAKq, 0.6iCo, 
O.Si^o, Ko, 2Ko, 3Ko, AKq, 5Ko, and lOKo (the lower 
curve) . 




T/K„ 



The chiral condensate S{T) (normalised to the that at 
T = 0) as a function of T (measured in the units of the 
Ko). 



which also demonstrates the same pattern of chiral sym- 
metry restoration. 



4. DISCUSSION 

In this paper we studied the GNJL model at finite 
temperatures. The effect of pair creation affected by 
both quark selfinteractions and temperature, is conve- 
niently described in the formalism of the Bogoliubov- 
Valatin transformations, parametrised by the order pa- 
rameter: the thermal chiral angle tpp{T). At T = 
spontaneous breaking of chiral symmetry happens as a 
result of such ^Pq pairs condensation in the vacuum. 
The resulting true (BCS) vacuum possesses a lower vac- 
uum energy than the trivial unbroken vacuum. At finite 
temperatures, two processes proceed in opposite direc- 
tions: confining interaction tends to condense more ^Pq 
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quark-antiquark pairs in the vacuum, that is to increase 
the chiral angle, while the effect of the temperature is 
a suppression of the latter. The balance of these two 
tendencies is encoded in the thermal mass-gap equa- 
tion (I13|) . Notice that, since the mass-gap equation is 
a nonlinear equation, then the interplay of these two 
effects is essentially nonlinear. For low T's the full ther- 
mal mass-gap equation can be split, however, into two 
parts: the nonlinear mass-gap equation which describes 
chiral symmetry breaking and the linear thermal equa- 
tion which describes its partial restoration due to the 
temperature. The resulting chiral angle is therefore just 
a sum of the zero-T angle ipp(T — 0) and the thermal 
correction dipp{T). Let us stress once more that such a 
disentanglement is only possible at T ^ Kq, where Kq is 
the natural energy scale introduced to the theory by the 
confining interaction. For high temperatures (T ^ Kq), 
chiral symmetry restoration is expected to take place, 
so the thermal mass-gap equation can be simplified un- 
der the assumption ipp{T) ^ 1. However, for T ~ Kq, 
the full nonlinear thermal mass-gap equation has to be 
considered and solved. 

As soon as the thermal mass-gap equation is solved 
and the chiral angle is found, it can be used in Bethe- 
Salpeter equations for hadrons or, alternatively, a sec- 
ond, generalised Bogoliubov-like transformation can be 
performed in order to diagonalise, for the hadronic sec- 
tor, the Hamiltonian of the theory therefore enabling 
us to build bound-state equations for mesons [5]. The 
proposed real-time approach is general and can be ap- 
plied to any microscopic model for confined systems. 
It allows one to investigate microscopically the proper 
balance between chiral symmetry breaking due to con- 
finement and chiral symmetry restoration due to the 
temperature. 

A brief comment on the temperature behaviour of 
the model is in order. An intrinsic feature of the GNJL 
model, with confinement provided by an infinitely rising 
potential, is that chiral symmetry breaking happens for 
any temperature T, thus the chiral symmetry restora- 
tion process being asymptotical. This is related to the 
fact that temperature effectively reduces the confining 
interaction strength, while even quite weak confinement 
is able to break chiral symmetry. Notice however that 
any realistic interquark interaction is expected to be 
temperature-dependent. In particular, lattice calcula- 
tions [TT] support the conjecture [H] that, while colour- 
magnetic vacuum fields do not change across the decon- 
finement phase transition, the QCD vacuum loses its 
confining colour-electric part. If this feature is incor- 
porated into the quark kernel, the resulting mass-gap 
equation will describe chiral phase transition. 



The small-T behaviour of the model also deserves 
special attention. From general considerations (see 
Ref. [13] for a review) one expects the following be- 
haviour of the chiral condensate at small T's: 

I](T)/I](0) = 1 - TV(8/2) + . . . . (23) 

However, since /tt ~ V^c, such corrections 0{T^) ap- 
pear only when the theory is considered in the order 
1/Nc or higher, that is beyond the mean- field approx- 
imation — see, for example, Ref. [14]. In the GNJL 
model the mean-field approximation corresponds to the 
BCS level, considered in this paper, while proceeding 
beyond the BCS approximation corresponds to the full 
diagonalisation of the Hamiltonian ([T]), including the 
: H4 : part — see Eq. ([3]) . The details of this diagonali- 
sation procedure can be found in Ref. [8] . The correction 
to the chiral condensate at small temperatures found in 
this paper is linear in T which differs from the behaviour 
obtained in other models (see, for example, the discus- 
sion of the NJL model in Ref. [M]). This should not 
come as a surprise however since the BCS vacuum of 
the model behaves like an infinite medium which is dis- 
turbed by the temperature, so that the response of this 
medium, in the leading order, is linear in the perturba- 
tion, that is in T. The medium has to possess rather 
peculiar features in order to respond quadratically or in 
higher powers of the temperature, while we do not ob- 
serve such properties of the BCS vacuum. However, be- 
yond the BCS approximation, the GNJL model is sub- 
ject to a severe realignment, so that it finally acquires a 
description entirely in terms of mesonic states [8j, with 
the massless pions being the lowest states in the spec- 
trum. Then, with the free energy calculated in terms 
of a pion gas, the standard approach can be applied 
to recover the usual behaviour ([23]) . Given a compli- 
cated form of the Bogoliubov-like transformation which 
diagonalises the Hamiltonian ^ up to the order 1/Nc 
(to be confronted to the 1/^/Nc diagonalisation at the 
BCS level), the cancellation of the linear in the temper- 
ature term in the chiral condensate does not look un- 
natural. Indeed, in the mesonic vacuum, only properly 
correlated pairs of mesons can be created or annihilated 
which appears to be a more involved process than ^Pq 
quark-antiquark pairs creation/annihilation in the BCS 
vacuum. Thus the mesonic vacuum is expected to re- 
spond weaker to the temperature increase, as compared 
to the quark vacuum. 

Notice that the approach suggested in this paper 
is very different from other, commonly adopted ap- 
proaches to investigation of the thermal properties of 
QCD. Indeed, while the "standard" procedure amounts 
to indirect tests of the vacuum through the properties 
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of the hadrons created on top of the vacuum, in our ap- 
proach the vacuum is probed directly. This makes im- 
portant further investigations of the interplay between 
our real-time approach and the standard imaginary- 
time approach to finite temperatures in confined sys- 
tems. We leave this problem for future publications. 
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